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Project  main  idea 

The  project  is  devoted  to  the  basic  research  and  establishes  possible  optical  ways  to  control  the  surface  plasmon 
polariton  in  graphene  layer.  A  system  comprises  the  graphene  ribbons  grating  placed  between  a  nematic  liquid  crystal 
(LC)  and  an  isotropic  dielectric.  Incident  light  wave  excites  a  plasmon  in  the  graphene  ribbons,  which  influences  the 
light  propagation  in  the  system.  The  grating  structure  of  the  graphene  monolayer  is  necessary  to  provide  the  plasmon 
excitation  because  there  is  a  huge  wave  vector  mismatch  between  the  graphene  plasmonic  wave  and  the  incident 
electromagnetic  wave.  In  graphene  structures  the  resonant  plasmon  frequency  depends  on  the  dielectric  properties  of 
layers  placed  above  and  below  the  graphene  sheet  or  ribbon.  In  our  case,  it  means  that  plasmon  frequency  must 
depend  on  the  LC  dielectric  permittivity.  The  LC  dielectric  permittivity  depends  on  the  LC  director  orientational 
state,  which  can  be  controlled  by  external  field.  Therefore,  a  propagation  of  the  surface  plasmon  polariton  in 
graphene  layer  can  be  tuned  by  reorienting  the  LC  director  using  external  electric  or  magnetic  fields. 

After  establishing  the  plasmon  freguency  dependence  on  the  LC  director  orientational  state  we  shall  study  a 
light  wave  propagation  across  the  system  under  consideration.  We  shall  calculate  the  absorption,  reflection,  and 
transmission  coefficients  and  show  that  control  of  the  orientational  state  of  the  LC  layer  enables  us  to  manipulate  by 
magnitude  of  the  absorption/reflection  maxima. 

These  results  will  allow  suggesting  a  new  type  of  the  graphene  micro-ribbon  grating  structures  to  control 
light  beams  propagation. 

Technical  approach 

Dependence  of  plasmon  frequency  on  the  LC  director  orientational  state. 

Schematic  of  the  hybrid  system,  LC  -  graphene  monolayer  -  isotropic  dielectric,  is  presented  in  Fig.  1. 


Fig.  1.  Schematic  of  the  hybrid  system:  LC  -  graphene  monolayer  -  isotropic  dielectric. 

We  use  the  following  designations:  £  1  is  the  LC  dielectric  tensor,  £2  is  the  isotropic  dielectric  permittivity,  G  is  the 
graphene  surface  conductivity,  y/,(p  are  the  LC  director  angles. 

We  write  the  Maxwell  equations  for  system  shown  in  Fig.  1  and  seek  their  solution  in  the  following  form: 

1)  for  the  wave  in  the  LC  (  z  >  0  ) 


E^(eA,+ezEtl)e^-m)-Q'\  H, 


2)  for  the  wave  in  the  isotropic  dielectric  (  z  <  0  ) 


1  U  2 


i(qx-cot)+Q2z 


(1) 


(2) 
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Substituting  (2),  (3)  into  Maxwell's  equations  we  take  into  account  that  the  LC  dielectric  tensor  has  a  general  form 


£xx 

£xz 

II 

£*y 

£yy 

£yz 

1 

V  £xz 

£yz 

£zz  y 

where  on  optical  frequencies 


sij=n02dij+[ne2-n02)didj. 


(3) 


(4) 


Here  dt  is  the  director  Cartesian  component,  no ,  ne  are  the  refractive  indices  for  ordinary  and  extraordinary  waves 
in  the  LC,  respectively. 

To  the  expressions  obtained  after  substitution  (1),  (2)  into  Maxwell's  equations  we  must  add  the  boundary 
conditions, 


E1  =  E] 

2x  Z=Q  lx  Z=Q 


•  K-M,.  <5> 


Demanding  the  boundary  conditions  (5)  to  be  satisfied  we  obtain  after  some  algebraic  transformations  a 
dispersion  equation  for  plasmon  wave  propagating  in  the  graphene 


v 


+  - 


I  2  C02Szz  2  £2® 

q - ^  A  q - ^2— 

c  V  c 


10 


cos. 


(6) 


were 


=  nQ2  +  (n/  -  n2 )  cos2  q>  sin2  y/, 
£zz  =n2  +[n2  -n2)  cos2  y/, 
sxz  =  [n2  -  n2 )  cos  (p  sin  y/  cos  y/ 


(7) 


and  the  LC  director  angles  are  shown  in  Fig.  1. 

The  graphene  surface  conductivity  G  in  (6),  in  general  case,  can  be  expressed  as  a  sum  of  two  terms, 
G int  er ,  where  the  first  term  corresponds  to  the  intraband  electron  transitions  and  the  second  term  gives 


G  =  G- 


int  ra  int  er 


contribution  from  the  interband  transitions.  For  fiCO  »  kBT  and  for  the  Fermi  energy  Ef  »  kBT  it  can  be 


approximately  written  as  [1], 


ie2 1  Er 


le 


-In 


2  EF\-h(o)  +  il rf 
2  EF  +ti((0  +  i / r) 


intra  j-2  /  ,  •  /  \  5  int  er  A  +. 

7rn  yco  +  i/T J  4 7rn 

In  the  nonretarded  limit  q  »  COS  /  C  the  dispersion  equation  (6)  simplifies  and  we  arrive  at  expression 


(8) 


ICOS ( 


q  =  - 


0  /  /  p  F  _  o2  F 

^xx^zz  bxz  '  ^2 


)• 


(9) 


It  is  seen  that  influence  of  the  LC  orientational  state  on  the  dispersion  equation  is  described  by  term 
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B(P>V')=yfe&r-£L  =  nen0y[si 


sin2  (p  cos2  y  - 


cos  <p- 


-  ( n0  /  ne  )2  sin2  <p  sin2  y 


(10) 


Reflection  and  transmission  coefficients. 

Consider  a  graphene  micro-ribbon  grating  in  the  xy-plane  placed  between  a  nematic  LC  layer  (anisotropic  top 
substrate)  and  an  isotropic  dielectric  substrate.  Each  micro-ribbon  in  the  grating  is  a  single  layer  graphene  with  the 
ribbon  along  the  y-axis,  where  d  is  the  micro-ribbon  width  and  A  is  the  grating  spacing  (Fig.  2).  A  plane 
monochromatic  wave  propagates  along  the  z-axis  from  the  side  of  the  LC  and  excites  the  plasmons  in  the  graphene 
micro-ribbons.  We  assume  semi-infinite  substrates  that  allows  us  to  neglect  the  effects  of  multiple  reflections;  in 
the  y  -  direction  the  system  is  infinite. 

To  simplify  calculations  we  suppose  the  LC  director  to  be  only  reoriented  in  the  xz-plane,  where  the  angle 
\f/  describes  the  director  deviation  from  the  z-axis.  We  also  set  the  magnetic  vector  of  the  incident  wave  to  be 
perpendicular  to  the  xz-plane  (TM-wave). 


Fig.2.  Schematic  of  the  graphene  micro-ribbon  grating  structure,  d  is  the  ribbon  width,  A  is  the  grating 
spacing,  y/  is  the  director  angle  with  the  z-axis,  8X  is  the  LC  dielectric  tensor,  s2  is  the  isotropic  substrate  dielectric 

permittivity,  <J  is  the  graphene  conductivity,  en  is  a  normal  to  the  graphene  plane,  and  k{  is  the  wave  vector  of  an 
incident  plane  monochromatic  wave. 

Because  of  the  strong  confinement  of  the  surface  plasmon-polaritons  (SPPs)  in  the  ribbons,  only  the  thin 
layer  of  the  LC  substrate  near  the  graphene  influences  the  SPPs.  Therefore,  we  can  set  the  LC  director  orientation  in 
the  whole  substrate  to  be  homogeneous  and  equal  to  the  director  orientation  of  this  layer  near  the  graphene.  The  LC 

optical  dielectric  tensor  can  be  written  in  the  form  £Uj  =  +  (if  —  ff  )n.n/ ,  i  =  x,  y,  z  ,  where  ni  denotes  the 

components  of  the  directors  =  (sin  ^/yOjCOS^)  ;  no  and  ne  are  the  refractive  indices  of  the  ordinary  and 
extraordinary  waves,  respectively  [2]. 

We  suppose  that  the  TM-wave  is  normally  incident  in  the  LC  on  the  graphene  grating.  As  it  follows  from 
the  Maxwell  equations,  electric  and  magnetic  vectors  of  this  wave  in  the  LC  take  the  form: 


E  =  (E  e  +  E  e  )  e 

i  V  ix  x  iz  z  ) 


-i(kt  z+co  t) 


H 


k; 


o 


E  -i(kiZ+*t)  E  ✓  I  )£ 

ix  y^iz  ^  Ixz  Izz  /  ix  ’ 


(ii) 


with  the  dispersion  equation 


(12) 
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For  reflected  and  transmitted  waves,  we  use  the  Fourier-Floquet  expansion  with  respect  to  the  coordinate  x. 
Satisfying  the  Maxwell  equations,  the  electric  and  magnetic  vectors  of  these  waves  can  be  written  as  follows: 

i)  for  the  reflected  wave 


E  = 


C  l  CO  £ 


lzz 


V  £\xx  S\xz  /  8lzz 


-  yT( k  +  £  k  )  e  —(s,  k  +  £  k  )  e  1  a  el("kn 

/  j  |_V  lzz  rn  1  xz  n  J  x  \  Ixz  m  Ixx  n  )  z  J  n 


x+krn  z-cot ) 


Hr  =  ^0  *x~SlJls\a  Z  an 

n 

with  the  dispersion  equation 


( knx+km  z-co  t ) 


ii)  for  the  transmitted  wave 


ev> 


k  = 


E  =- 


2A 


8  Ixx 


'Ixz 


'lzz  J 


( co2  k2^ 


\c  J 


"1X2 


'lzz 


CO £2  n 

H  t=S0C^tY,bneHK 


x—ktn  z~co  t ) 


X-kfn  Z-COt) 


(13) 


(14) 


(15) 


with  the  dispersion  equation 


k  -  \0_e  _k2 
Ktn  ~  H  ^2  ^2  Kn 


(16) 


Here  k  n=  2 7rn  /  A  ,  n  is  the  number  of  a  Floquet  spatial  harmonic. 

We  define  the  current  density  in  the  plane  of  the  graphene  grating  as js  =<j(x)Etx( z  =  0)ex, 
where  <j(x)  =  o  within  each  micro-ribbon  and  cr(x)  =  0  in  the  gaps  between  micro-ribbons. 
Substituting  Etx(z  =  0)  from  eq.  (15),  we  obtain 


MO- 


CD 


5X 


b  0e 


(17) 


To  obtain  equations  for  coefficients  an ,  bn  of  the  Fourier-Floquet  expansions  (13),  (15)  we  use  the 
boundary  conditions: 


en  x(H,.  +Hr  -Hr  )-/,]| z=0  -ex  =  0, 

z  -1  (to) 

_en  X  (E,-  +  Er  -  Ef  )]  I  z=0  -ty  =  0. 

Now  we  substitute  eqs.  (11),  (13),  (15)  for  electric  and  magnetic  vectors  and  eq.(17)  for  the  current 
density  js  into  eqs.  (18).  Then,  using  the  dispersion  equations  (12),  1(4),  (16)  we  arrive  at 
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-2  in  mx 

e  x 

Multiplying  eqs.  (19),  (20)  by  ,  where  m  is  an  integer  and  integrating  these  equations  over  x  in 

the  range  [0,  A]  we  obtain  the  set  of  equations  for  coefficients  an,  bn: 


In  the  infrared  and  terahertz  spectral  regions  the  intraband  contribution  to  the  graphene  conductivity 
dominates.  Then,  in  the  random  phase  approximation  the  2D  conductivity  of  graphene  can  be  written  as  [3] 


2  77 

e  E¥  ti 
xh2  cor  +  i  ’ 


(23) 


where  EF  is  the  Fermi  energy  of  graphene,  T  is  the  electron  relaxation  time,  e  is  the  electron  charge. 

The  reflection  and  transmission  coefficients  are  defined  as 

R  =  |Re(Er  xH*)|/|Re(E;  xH*)| ,  T  =  |Re(Ef  xH*)|/|Re(E;  xH*)| .  (24) 
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The  graphene  absorption  coefficient  isA  =  l  —  (R  +  T)  .  Numerically  solving  equations  (21)-(23)  the  coefficients 

an,bndiVQ  obtained.  Using  formulas  (13),  (15)  to  find  the  electric  and  magnetic  vectors  of  the  reflected  and 

transmitted  waves  and  after  their  substitution  into  formula  (24),  we  calculate  the  reflection,  transmission,  and 
absorption  coefficients. 


Technical  progress  overview 

In  Fig.  3  we  show  dependence  of  quantity  E^(p,y/^  on  the  LC  director  angles  (p,\f/  for  two  pairs  of  the  refractive 
indices  of  the  ordinary  and  extraordinary  waves  in  the  LC,  no ,  fle  . 


Fig.  3.  E{^cp,y/)  versus  (p,y/  :  (a)  n0  =1.53,  ne  =  1.94 ;  (b)  n0  =1.94,  ne  =  1.53  . 
Below  in  Fig.  4  we  also  show  for  illustration  the  cross-sections  of  function  E  {(p,  l//^j  by  vertical  planes. 
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Analyzing  formula  (10)  (or  Figs.  3,  4)  we  can  conclude  that  ratio  i?max  /  Bm-n  responsible  for  interval  of  tuning  the 
plasmon  wave  vector  (or  frequency)  by  the  LC  director  reorientation  equals  ne  /  no  if  Tle  >  no  or  no  /  Tle  if 
fl0  >  fle  .  It  is  necessary  to  note  that  in  the  retarded  region  ( q  <  COS /  c)  where  instead  of  formula  (9)  we  have  to 
use  formula  (6),  the  small  correction  to  this  estimation  can  arise  from  Szz  in  denominator  of  formula  (6). 

As  material  of  the  nematic  layer,  we  use  the  high  birefringence  (An  «  0.41)  LC  mixture  W1791 
with  n0  ~  1.53  ,  ne  ~  1.94  at  X  —  1.064  pm  [4].  The  birefringence  of  W1791,  which  is  known  in  visible  and 

near-IR  regions  of  spectrum,  was  extrapolated  to  mid-IR  and  THz  frequencies  using  the  extended  Cauchy  dispersion 
formulas  obtained  for  the  LC  mixture  E7  [5].  In  our  study,  we  chose  either  a  low  dielectric  constant  material  (e.g. 
hexagonal  boron  nitride,  h-BN,  S2  =  Sh_BN  =  3 )  [6,  7]  or  a  high  dielectric  constant  material  (silicon,  Si, 

8 2  ~  8 si  =  1 1  -  7  )  [8]  as  the  isotropic  substrate  material. 

The  physical  parameters  of  the  graphene  grating  in  our  simulations  are  as  follows:  the  grating  spacing  was 
fixed  at  A  =  1  pm,  the  micro-ribbon  width  was  varied  as  =  0.3  A  ,  0.5A  and  0.7  A  .  For  evaluation  of  the 

electron  relaxation  time  we  used  the  formula  T  =  jdEF  /(evF)  ,  where  jd  is  the  carrier  mobility,  vF  =  3  x  106m-s_1 
was  the  Fermi  velocity  in  graphene  [10].  Setting^  =  0.64  eV  [9]  and  the  carrier  mobility  jd  —  0.5  m2/(V-s)  we 

get  the  carrier  scattering  time  T  —  0.32  ps.  The  simulations  were  performed  in  the  wavelength  range  14.5-55  pm 
(20-5.45  THz)  with  a  temperature  T  =  300  K  . 

The  convergence  of  the  computational  procedure  used  with  respect  to  higher  Fouquet  harmonics  is 
necessary.  In  order  to  ensure  the  required  accuracy  of  calculations,  it  was  sufficient  to  select  the  number  of 
harmonics  N  >  400  in  eqs.  (21),  (22).  Results  of  our  numerical  calculations  of  absorption,  reflection,  and 
transmission  spectra  of  the  graphene  grating  are  shown  in  Figs.  5  and  6. 

Fig.  5  illustrates  the  change  of  the  absorption,  reflection,  and  transmission  spectra  of  the  graphene  grating 

for  two  limiting  values  of  the  LC  director  angle  y/  —  0  and  iff  =  90° ,  and  a  ribbon  width  to  grating  spacing 
ratio  d  /  A  =  0.5  .  Curves  1  and  2  correspond  to  the  cases  of  an  isotropic  substrate  with  either  a  low  or  a  high 
dielectric  constant,  respectively.  Maxima  in  the  absorption  and  reflection  spectra  are  related  to  the  excitation  of  the 
plasmons  in  the  graphene  micro-ribbons.  Within  our  frequency  range  two  plasmon  peaks  in  curves  2  are  observed. 
The  ratio  of  the  resonant  frequencies  corresponding  to  these  peaks  is  equal  to  2,  which  agrees  with  the  results 
obtained  in  paper  [9]  for  the  graphene  grating  placed  between  two  isotropic  dielectric  substrates. 

As  we  can  see  from  Fig.  5,  values  of  these  maxima  depend  on  the  LC  director  orientation.  In  particular,  the 

rotation  of  the  LC  director  by  90°  leads  to  a  change  of  the  absorption  maximum  value  by  approximately  14%  when 
£*2=3  (h-BN  substrate)  and  by  16%  when  £*2=11.7  (Si  substrate). 


Fig.  5.  Absorption  (a),  reflection  and  transmission  (b)  spectra  of  the  graphene  micro-ribbon  grating  at  different 
angles  of  the  nematic  director:  y/  —  0°  -  solid  line,  y/  =  90°  -  dash-dot  line.  Substrate  dielectric  constant: 

£*2  =  3  -  curves  1,  £*2  =  1 1 .7  -  curves  2.  Ratio  d  /  A  =  0.5  . 
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In  Fig.  6,  the  change  of  the  absorption,  reflection,  and  transmission  spectra  of  the  graphene  grating  at  the 
two  limiting  values  of  the  LC  director  angle  y/  is  shown  for  two  values  of  the  ribbon  width  to  grating  spacing  ratio, 

d  /  A  =  0.3  and  d  /  A  =  0.7  .  In  this  case,  the  rotation  of  the  LC  director  by  90°  leads  to  the  change  in  the 
absorption  by  approximately  14%  when  d  /  A  =  0.3  and  by  18%  when  d  /  A  =  0.7  .A  comparison  of  the  results 
presented  in  Figs.  5,  6  shows  that  the  effect  of  the  LC  director  rotation  increases  with  an  increase  of  the  isotropic 
substrate  dielectric  constant  and  the  ribbon  width  to  grating  spacing  ratio  d  /  A  . 


Fig.  6.  Absorption  (a),  reflection  and  transmission  (b)  spectra  of  the  graphene  micro-ribbon  grating  at 
different  angles  of  the  nematic  director:  y/  =  0°  -  solid  line,  \f/  =  90°  -  dash-dot  line.  Ratio  d  /  A  =  0.3  -  curves 

1,  d  /  A  =  0.7  -  curves  2.  Substrate  dielectric  constant  s2  =11.7. 

On  scientific  results  of  the  project  one  paper  is  published,  one  paper  is  prepared  for  publication,  and  four 
reports  on  international  conferences  were  delivered. 
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Current  status  of  the  project. 

Current  technical  status  of  the  project:  on  schedule. 

Summary  of  personnel  commitment. 

Professors  of  Taras  Shevchenko  National  University  of  Kyiv  V.  Yu.  Reshetnyak,  I.  P.  Pinkevych  and  Dr.  V.  I. 
Zadorozhnii  fulfilled  the  tasks  during  the  reported  period. 

Prof.  V.  Yu.  Reshetnyak  -  developed  the  mathematical  model  of  the  hybrid  system  with  the  LC  layer, 

-  obtained  expressions  for  plasmon  frequency  in  the  hybrid  system  with  the  LC  layer, 

-  calculated  the  graphene  plasmon  frequency  in  dependence  on  the  orientational  state  of  the  LC  layer, 

-  obtained  equations  describing  light  propagation  in  the  hybrid  system  with  plasmon  excitation  in  the  graphene, 

-  developed  numerical  methods  for  calculation  coefficients  of  light  reflection  and  transmission  in  the  hybrid  system, 

-  optimized  the  hybrid  system  parameters  for  active  control  of  light  propagation  and  reflection. 

Prof.  I.P.  Pinkevych  -  obtained  and  solved  equations  for  the  director  and  electric  field  in  the  LC  cell, 

-  studied  the  influence  of  anchoring  on  the  LC  cell  boundaries  on  the  director  in  the  cell  bulk, 

-  studied  the  control  of  plasmon  frequency  by  DC  field  and  incident  light  field, 

-  obtained  expressions  for  coefficients  of  light  reflection  and  transmission  in  the  hybrid  system, 

-  studied  influence  of  the  LC  layer  orientational  state  on  light  reflection  and  transmission  in  the  hybrid  system. 

Dr.  V.I.  Zadorozhnii  -  fulfilled  numerical  calculations  of  light  reflection  and  transmission  coefficients  in  dependence 
on  the  anchoring  conditions  on  the  LC  cell  boundaries. 

Variations  in  the  scheduled  amounts  of  efforts. 

Prof.  V.  Yu.  Reshetnyak  and  Prof.  I.P.  Pinkevych  fulfilled  the  following  tasks  planned  initially  for  Dr.  V.I. 
Zadorozhnii: 

calculation  of  the  graphene  plasmon  frequency  in  dependence  on  the  orientational  state  of  the  LC  layer 
(Reshetnyak), 

develop  numerical  methods  for  calculation  coefficients  of  light  reflection  and  transmission  in  the  hybrid  system 
(Reshetnyak) 

study  the  influence  of  anchoring  on  the  LC  cell  boundaries  on  the  director  in  the  cell  bulk  (Pinkevych). 

Description  of  travels. 

Travel  of  prof.  V.  Yu.  Reshetnyak  for  participation  in  Conference  on  Photorefractive  Materials,  Villars,  Switzerland, 
16-19  June  2015.  Report  was  delivered  on  scientific  results  of  the  project. 

Travel  of  prof.  V.  Yu.  Reshetnyak  for  participation  in  9th  International  Congress  on  Advanced  Electromagnetic 
Materials  in  Microwaves  and  Optics,  Oxford,  UK,  7-12  September  2015; 

Travels  of  prof.  V.  Yu.  Reshetnyak  and  prof.  I.  P.  Pinkevych  for  participation  in  16th  Topical  Meeting  on  the  Optics 
of  Liquid  Crystals,  Sopot,  Poland,  13-18  September,  2015.  Report  was  delivered  on  scientific  results  of  the  project. 
Travels  of  prof.  V.  Yu.  Reshetnyak  and  prof.  I.P.  Pinkevych  for  participation  in  International  conference  "Flowing 
Matter  2016",  Porto,  Portugal,  11-15  January  2016.  Report  was  delivered  on  scientific  results  of  the  project. 

Information  about  major  equipment  and  materials  acquired,  other  direct  costs,  related  to  the  project. 

There  was  not  purchase  of  equipment  or  materials,  and  other  direct  costs  have  not  been  spent. 


Table  Redirection 


Reference  New  requested  category,  Requested 
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(1)  (2)  (3) 
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1.1  Grant  -  FWS  |  777g  ?5 
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